Вестник ДГТУ, 2010. Т. 10. №4(47) 








ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ 


УДК 539.3 
Д.А. ПОЖАРСКИЙ, А.А. МОЛЧАНОВ 


АСИМПТОТИЧЕСКИЕ РЕШЕНИЯ СМЕШАННЫХ ЗАДАЧ 
ДЛЯ УПРУГОЙ ПОЛОСЫ И КЛИНА 


Получены асимптотические решения двух новых задач теории упругости со смешанными граничными усло- 
виями. Первая задача — о трещине сдвига в полосе при разных вариантах граничных условий на внешних 
гранях полосы. Вторая — контактная задача для трехмерного клина при действии сосредоточенной силы вне 
области контакта на ребре клина и скользящей заделке другой грани. 

Ключевые слова: теория упругости, смешанные граничные условия, асимптотики. 


Введение. Актуальность задач механики разрушения [1] не вызывает сомнений. Однако до сих 
пор были исследованы только задачи о трещине нормального отрыва в упругой полосе [2]. В на- 
стоящей работе строится асимптотическое решение задачи о трещине продольного сдвига в по- 
лосе. Задачи о трещинах тесно связаны с контактными задачами (в обоих типах задач возникают 
смешанные граничные условия) [3—5]. Реальная поверхность часто имеет микронеровности, по- 
этому контакт может происходить одновременно на нескольких участках. Для изучения дискрет- 
ного контакта (взаимодействия нескольких штампов) важно заранее получить решение задачи 
о действии дополнительной сосредоточенной силы вне области контакта упругого клина [4, 5]. 

1. Трещина сдвига в полосе. Постановка задачи. В декартовых координатах х, у рассмотрим 
упругую полосу {хе (—,0); уе [0,21] } толщины 2й. Предполагаем, что в рамках линейной теории 
упругости материал полосы характеризуется двумя параметрами упругости: С — модуль сдвига 
иу- коэффициент Пуассона). Пусть обе грани полосы у = 0, у = 2й жестко заделаны (вариант А), 
находятся в условиях скользящей заделки (вариант Б) или свободны от напряжений (вариант В). 
В середине полосы при у = й имеется трещина продольного сдвига (разрез), занимающая область 
ХЕ (-а,а), нагруженная по обоим берегам заданной касательной нагрузкой т(х). Полагаем, что 
нормальное напряжение оу = 0 при хе (—0,0), у = й [1]. Требуется определить скачок касатель- 
ного перемещения в области трещины и(х) = +и,(х,й-0). Затем может быть определен коэффици- 
ент интенсивности касательного напряжения в кончике трещины, ответственный (согласно крите- 
рию разрушения) за ее дальнейшее распространение. 

Асимптотическое решение. В силу одинаковости напряженно-деформированного состояния 
достаточно рассматривать лишь половину полосы 0<у<й. Для вывода интегрального уравнения 
относительно функции и(х) рассмотрим вспомогательные краевые задачи с граничными условия- 
ми вида 


у=0: А) и, =0, и, =0; (1) 
у=0: Б) т,, =0, и, =0; (2) 
у=0: В) т,,=0, в, =0; (3) 
у=й: &,=0; и, =—и(х). (4) 


Здесь и(х) - временно известная функция, равная нулю при хе (-—а,а). Для решения краевых 
задач (1)-(4) воспользуемся общим решением дифференциальных уравнений равновесия Ламе 
плоской задачи теории упругости в форме Папковича-Нейбера и преобразованием Фурье. Затем 
для получения интегрального уравнения в области трещины используем граничное условие 
(с заданной правой частью) 


у (х, п) = —(х), (5) 
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при помощи которого относительно функции и(х) выводим интегральное уравнение (для простоты 
ограничимся случаем х(х) = т = соп$®: 





орел (ка): (6) 
С С 

(5) = | ИК) сов в, 0=—; (7) 
]1-у 


0 
где для варианта А 


ксв21 + 21? +0,51+к?). 





ко к $12 + 21 | (8) 
а для вариантов Б и В соответственно 
2 2 
р. И" 1—1 ) (9) 
сВ2#+1 $621 -— 21 
Решение уравнения (6) подчиним условиям 
и(+а) = 0. (10) 


Отметим сходство уравнения (6) с уравнением (8.2) [2] для задачи о трещине нормально- 
го отрыва в полосе. Это позволяет применить для решения уравнения (6) регулярный асимптоти- 
ческий метод эффективный для относительно толстой полосы [2]. Интегрируя по частям в (6) при 
учете (10) и вводя безразмерные величины 








х =; и=-; Л(х )=-_ибх); т =; (11) 
а а ах 0 
получим уравнение (звездочки далее опускаем) 
1 
—х 
| лещ - } =-лье (ху: (12) 
—1 
Ко (5) = | Г) ут 5 @. (13) 
0 
Выделим из ядра (12) главную часть в форме 
о | [2(0) — Изв 6 4. (14) 
5 
0 
Гладкая часть ядра (5) при |5|<2 разложима в абсолютно сходящийся ряд 
Н()= У е,5""Н; (15) 
п=0 
ра [о ПИ"! @ (п=0,1,2,,..). (16) 
(21 +1)! 5 


Несколько первых значений постоянных (16) для вариантов А, Б и В даны в табл.1 
(при т = 0,3). 

















Таблица 1 
Значения постоянных (16) 
Задача @ е, [=> 6 
А 0,9429 —0,5435 0,2930 —0,1301 
Б 0,4112 —0,3551 0,1540 —0,05448 
В —1,3337 0,9667 —0,5153 0,2286 
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Регулярное асимптотическое решение уравнения (12) имеет вид [2] 


(о) = -- РФ); (17) 


М -х2 


2 
Ф(х, и) =1 0 а — " 1 + 
а 4 8 2 


у еб + Зеоел , 15е2 , Зее,  5е> №, 5е> 4 1 ВО 1 (18) 
8 16 8 4 2 2 цб из 


Из (17), (18) найдем отнесенное к а касательное раскрытие трещины: 


х 2 
и = [| убощко х? | “о. -[ Ав - - + 
а 
—1 


22 |4 8 2 


3 
О Пере, 65е› + 2061 25е› ый 5е> 4 1 ие 1 (19) 
8 16 72 4 18 12 цб и 


Размерный коэффициент интенсивности касательных напряжений (КИН) в кончике трещи- 
ны, отвечающий за дальнейшее распространение трещины, имеет вид (величина х - размерная) 




















Ки=- Ши 0./2п(а-х) т (20) 
х—>>а-0 Ах 
Безразмерный КИН на основании (19), (20) имеет вид (величина т - размерная) 
К | т 
а о. +| 90 9е1 - т 5ее, 25е> ы +0 г | (21) 
т/ап 24 4 8]и 8 16 8 |н н 


Отсюда, используя табл.1, получим для варианта: 











К 
И пы 0,4715 08337 1500 , а | } 5 
туал п п п п 
К 
ее мес 0.2056 ов 0.6269 й й .) бя 
туат п п п п 
К 
да и. О 0.6428 0,6982 ь а 1 } и 
туал п п п п 


Формулы (22)-(24) рекомендуется использовать при и > 2. При малых значениях д можно 
использовать сингулярный асимптотический метод [2]. 
2. Контактная задача для клина. Постановка задачи. В цилиндрических координатах и, ф, 2 
рассмотрим трехмерный упругий клин {те [0,5], фе [0,0], 2Е (—00,0)} угла раствора с, ось 2 на- 
правлена по ребру клина. Материал клина характеризуется модулем сдвига и коэффициентом Пу- 
ассона. Грань ф = 0 находится в условиях скользящей заделки. В грань ф = о без перекоса вдав- 
ливается эллиптический в плане штамп с плоским основанием, к которому приложена сила Р. Об- 
ласть контакта — эллипс © = {(7-а)”/с*+2/Ь*}, Ь>с, а>с. В начале координат г = 2 = 0 к грани 
ф = а приложена нормальная сосредоточенная сила О. Для простоты считаем задачу симметрич- 
ной по координате 2. Осадка штампа равна 6. При заданных величинах си, С, у, а, Ь, си д требует- 
ся определить контактное давление под штампом о. = -94(7,2), (7,2) © и затем найти силу Р. 
Асимптотическое решение. Сила О вызывает в заданной области контакта <> дополнительное 
нормальное перемещение [3] 


ио(т, 0,2) =-_ 94002), ве (25) 


210? + 22 1-у 
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© Е 
ры [* > И’(и)Е (и) = Ш ен] Чи; (26) 
п г 
0 
р 028 ое ЕЯ (27) 
2а +;щ 2а $62ам + изш 2а 
Функция Е(и) удовлетворяет интегральному уравнению Фредгольма второго рода (0<и<оо) 
Е(и)--2%) | Ки, у)Е (ау =п/2 (- 2У)Ки,0); (28) 
0 
Ци.) = сви в (У) [ ЕО (29) 
2 2 ы (сплЕ + сЬли)(сЬлЕ + слу) 
со 02 
2 
ЕО ли [и (Оф > ЗЕ сои зт” 24а (30) 
2 о свлё + сБли сВ204 — с0$ 44, 
В частном случае о=л/2 на основании формул (25) - (30) получим 
ии, /2, 2) =, (31) 


л0\/ "+ =? 


что соответствует двум одинаковым силам на упругом полупространстве [3]. 
Используя (25) и известную функцию Грина для трехмерного клина [3], получим относи- 
тельно 4(7,2) следующее интегральное уравнение: 


А 
т показ О о) | (,2) ей; (32) 
0 2 2 
[9] 2п0\и"- +2 
зы . р_ 2 2.1/2. 
Коро ИА) В =[(и-х) +=) ] > (33) 


Е(х, у,г,2)= = | | (звли(И/ (и) — свли)К (Вх) + 
00 


+ 256 — ИШЕ(и, Вх} К д, (Ви) соз В(2 — у)аВаи. (34) 


Здесь К„(г) - функция Бесселя, а функция Е(и,Вх) удовлетворяет интегральному уравне- 
нию Фредгольма второго рода (0<и<оо) 


Е(и,Вх) = - 25) [ Ци, УБЕ. вл) + св Ку (Вх) (35) 
0 


Отметим, что перемещение (25) получается из ядра интегрального уравнения (32) при 
х—>0, у = 0 (с точностью до множителя). При о=л/2 ядро интегрального уравнения (32) совпадает 
с ядром интегрального уравнения симметричной контактной задачи о двух штампах на полупро- 
странстве. 

Для получения асимптотического решения введем безразмерные величины (штрихи далее 
опускаем) 

г-а 2 а с б К 

= =—, == =—, А'=—, В'= ВБ; 36 
Ь г ь р В'=В (36) 
4'(",2')=9(+,2)/0, О'=О К2л0Б?), Р'=Р К0Б?), ©'<> ©. (37) 











) ) 
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Уравнение (32) в обозначениях (36), (37) принимает вид ( (",2)е О) 


ОА (и), 2) 


(+)? + =? 


Для решения интегрального уравнения (38) применим регулярный асимптотический 
метод [3], эффективный при больших значениях параметра /,, т.е. вдали от ребра клина. 

С использованием известных интегралов разложим гладкую часть ядра уравнения (38) в 
ряд по отрицательным степеням ^.. Имеем (/, —> о) 





| а(х, Я + РА, у г+А, 264 =2л] 5 (38) 
о 








А А Аз (х? +7) + А. (у-2) + А 
А (39) 
А, 7. А ры 
А 3 а+к 
А а о в бе 40 
1 =@0 +Ко, 45 5 3340 5 | (40) 
Ч 90 
Ад =—-к›, == -а чу, 41 
4=> ет 1 (41) 
где (и = 0,1, т = 0,1,2): 
© | в о 
а, =——— | [®лиЙ’ (и) -— 1(0,25 +и^) "аи; (42) 
р 
ежыо | ИЕ, (маи. (43) 
0 


Значения функций Р„(и’), где и, — узел квадратурной формулы Гаусса, определяются 


как значения при и = и, функций, удовлетворяющих интегральным уравнениям Фредгольма вто- 
рого рода (0<и<о) 





ея | сН(пу/2)„ (их, |. 
Е„(и)= (1-25) ] Ци, у) [А а | (44) 
19(и,у)=16, [5 (и,у)=(1+(и+ у)’ +(и-у)?); (45) 
[ (и, у)=(-—(и+у)?)И+(и- у)?)+(3—(и-у)?)а+(и+у). (46) 
Аналогично разложим функцию в правой части уравнения (32). С учетом того, что 
(^, >) 
2 2 

И +} (47) 

("+22 +22 ^ № 2% йо 





("+ ^)? + 22+ 2:2 
(+)? +22+|2| а (1 +0]; 





с05| и Ш —. 48 
й ИА 942 А, И (18) 
получим разложение (4, —> с): 
2 2 
Ар(’+А,2) _ В и ыы +0 + (49) 
(+)? +22 ^^ ^ ^ 
А+Ь +Ь 
В! =-В> = В =А+Ы, и. (50) 
р | И Ем) и?" аи (п=0)). (51) 
пу 2 
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и (49), рассчитанные при у = 0,3 и а = лА/б. 


Здесь постоянная А определяется формулой (27), а функция ЁР(и) удовлетворяет инте- 
гральному уравнению Фредгольма (28). 
Разложения (39) и (49), если их продолжить, равномерно сходятся в области (х,у),(7,2)Е © 
при ^, > тах(о ',1+е). 
В табл.2 приводятся значения независимых коэффициентов, входящих в разложения (39) 


























Таблица 2 
Значения независимых постоянных (40), (41), (50) 

К 2 3 4 5 6 7 8 
— А! 0,2867 —0,5000 —0,2701 0,1732 0,3183 0,2467 0,1502 
— —0,1400 —0,1250 —0,06667 0,05354 0,09284 0,07298 0,04440 

Аз —0,007466 | -—0,06250 —0,04353 0,007344 0,02653 0,01766 0,01036 
—/5 | -0,006701 | -0,2500 0,1367 0,06615 0,1326 0,1008 0,06138 

В 1,843 2,000 1,494 0,3694 0 0,1944 0,5160 
_ В 1,096 1,000 0,7944 0,1925 0 0,09939 0,2625 


























При учете (39) и (49) асимптотическое решение уравнения (32) при больших значениях /, 
ищем в виде 





3 
ан) +0. (5%). (52) 
п=0 ^. ^ 


Приходим к цепочке интегральных уравнений с ядром А ' для последовательного опреде- 
ления функций из разложения (52). Каждое из уравнений этой цепочки решается в явном виде. 
В результате получим 








___б _ __- О! в | и. 0, 
Е Е 5(х, у) = 1 =2 У, (53) 
_ бАи м > А! О зе? х . 
От = ОВ» 5 К (т = 1,2,3,4), 42(х, у) = => К? = ' 
1 1|_ РО, _ 8(43=* + 44) 2 ь 
Е | К? ЗК у Ра + ОзЕ "> (54) 


2 


А. Оле?х 4х? 4 
2 + (От - Оу зе’ -ои, е=\1-=? , 


К(Е-=?К) \ 
где К = К(е), Е = Е(е) - полные эллиптические интегралы; 














уи==(=е”)К” +20 -е”“)КЕ-ЗЕ”; (55) 
и, = (1-е?) (1,5 -е?)К? -2(1-е?)КЕ + 0,5Е?; (56) 
> = (1,5 — 0,5е? )К? -2КЕ +0,5Е?, уу =К-(1+0,5е?)Е; (57) 
М =\Ую =Е-@-0,5е?)К, у! =(-е?)? К -(-1,5е?)Е. (58) 
Для вдавливающей штамп силы получим выражение 
3 
Р 
Ра [аоьз)аху = 2». “+О : (\. >); (59) 
т х 
о п=0 
б О АО 
а (60) 
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а - А20, 6(4Е?+ 4.) 2 
Ри _ 


Е РЕ ЗК ри оу) 


4 
+” ау — Озу о) + ОЕ Чо — Оли]. (61) 


Расчеты по формулам (52)-(61) показывают, что при возрастании дополнительной силы 
О, приложенной вне штампа, снижается как контактное давление, так и сила, вдавливающая 
штамп. При возрастании значения /, (область контакта отдаляется от ребра) контактное давление 
и вдавливающая сила стремятся к известным значениям [3]. 
Выводы. Полученные формулы (22)-(24) для коэффициента интенсивности касательных напря- 
жений в кончике трещины важны для анализа возможности дальнейшего распространения тре- 
щины. Этот анализ показывает, что трещина будет наименее опасна для граничных условий ва- 
рианта А (жесткая заделка граней полосы) и наиболее опасна для варианта В (свободные грани 
полосы). Решение задачи о действии сосредоточенной силы вне области контакта на грани клина 
может быть использовано в дальнейшем для сведения контактной задачи о взаимодействии не- 
скольких штампов на грани клина к системе интегральных уравнений Фредгольма второго рода 
(метод Андрейкива - Панасюка) [5]. 

Работа для клина поддержана грантами РФФИ 08-01-00003, 09-01-00004. 
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АЗУМРТОТТС $ОГОТТОМ$ ОЕ МТХЕО РКОВЕЕМ$ РОК ЕТАЗТТС 5ТВТР АМО ММЕОСЕ 
АзутроНс 5оиНопе Гог {мю пем! {Пеогу еазНсКу ргоБМетз мА плхед Боипдагу сопаоп$ аге аейуеч. 
Тре Яг${ ргоБет [5 аБоцЁ а $Пеаг сгаск т а $ р Гог АаТегепе {уре$ оЁ Боипдагу сопаопз оп {Ле ощег 
$ р Гасез. ТНе зесопа опе {5 {Не сощасЕ ргоШет Гог а {Агее-Айтепзюпа! ммедде мтеп а сопсепгаеа 
Гогсе 15 асупа оцёае {Пе сопкасЕ дотат оп {Пе медде едде Гог $!4тд зиррой оп Не о{ег медде 


Гасе. 
Кеу мгог4$: Пеогу ог е!азНсКу, птхей Боупдагу сопа®оп$, азутрю{сс. 
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